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I .  INTRODUCTION 


This  report  documents  the  form-functions  developed  for  the  interior 
ballistics  computer  code  Interior  Ballistics  of  High-Velocity  Guns 
(IBHVG) ,  a  FORTRAN  program  used  in-" the  Interior  Ballistics  Division  of 
the  Ballistic  Research  Laboratory.  IBHVG,  based  on  the  Baer-Frankle 
code1,  will  itself  be  the  subject  of  a  forthcoming  report. 

A  "form-function"  is  any  computational  scheme  for  determining 
either  surface  area,  volume,  or  fraction  burned  for  a  burning  propellant 
grain  at  any  instant  from  ignition  to  extinction  based  on  initial  dimen¬ 
sions  and  depth  burned.  Hither  surface  area  or  fraction  burned  is 
generally  required  in  Interior  ballistic  calculations.  The  IBHVG  code 
does  not  use  surface  area  so  the  form- functions  described  below  omit 
area  derivation.  Also,  the  traditional  assumptions  are  made  that  all 
grain  surfaces  erode  at  the  same  rate  and  all  the  perforations  are  of 
equal  diameter. 

A  variety  of  approaches  to  a  number  of  grain  geometries  have  been 
taken  in  the  past2**.  Stals,  for  example,  not  only  deals  with  detailed 
analyses  of  different  grain  types,  but  contains  a  wealth  of  background 
information,  history  of  attacks  to  the  problem,  ard  an  extensive  biblio¬ 
graphy.  While  this  report  does  not  attempt  to  supplant  previous  efforts, 
it  presents  some  novel  treatments  which  may  prove  applicable  to  codes 
other  than  IBHVG.  For  example,  the  analyses  of  the  seven-  and  nineteen- 
perforated  cases  are  based  on  a  method  of  dealing  with  general  cylindri¬ 
cal  grains*.  This  method  enables  handling  grains  with  unequal  webs. 

The  FORTRAN  coding  embodying  the  methods  of  calculations  will  be 
discussed  in  parallel  with  the  development  of  all  algorithms  for  oase  of 


1 

P.  G.  Baer,  and  J.  M.  Frankie,  "The  Simulation  of  Interior  Ballistic 
Performance  of  Gum  by  Digital  Computer  Program USA  Ballistic  Research 
Laboratories  Report  1183 ,  USA  Ballistic  Research  Laboratories,  APG,  MD, 
December  1962 . 

2 

V.  Comer,  "Theory  of  the  Interior  Ballietios  of  Guns,"  John  Wiley  A 
Sons,  New  York ,  1060,  pp,  30-36, 

z 

J .  Stale ,  " Form-Functions  for  Multicomponent  Propellant  Charges 
Including  Inhibited  Grains  and  Sliver  Bum,"  Materials  Reosaroh  Labor¬ 
atories  Technical  Note  371,  Materials  Research  Laboratories,  Maribymang , 
Victoria,  Australia,  September  1976, 

^F.  R.  Lynn ,  "Development  of  General  Form-Functions  for  Multiperfor- 
ated  Cylindrical  Propellant  Grains Ballistic  Research  Laboratory 
Memorandum  Report  ARBRLMR03014,  Ballistic  Research  Laboratory,  USA 
ARRADCOM,  APG,  MD,  March  1980. 
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understanding.  Also,  FORTRAN  computational  hierarchy  is  followed  through¬ 
out  the  text.  A  comprehensive  listing  of  the  referenced  subroutines 
will  be  found  in  APPHNDIX  A. 


II.  SUBROUTING  ORGANIZATION 

The  result  of  the  form-function  calculations  in  IBHVG  is  the 
fractions  burned,  at  any  instant,  for  the  M  propellants  comprising  the 
charge,  where  1  <  H  <_  5.  At  each  time  step,  a  call  is  made  to  a  single 
subroutine,  FORMT,  which  either  performs  the  required  operations  or,  in 
the  case  of  seven-  or  nineteen-perforated  grains,  calls  auxiliary 
routines  GBNIS  and  GENOS.  APPENDIX  A  contains  listings  for  these  three 
subroutines.  v 

APPENDIX  B  notes  all  important  input,  output,  ana  control  variables 
used.  Most  are  passed  through  COMMON,  although  a  few  are  local.'  No 
other  COMMON  variables  pertain  to  FORMT  operation. 

Subroutine  FORMT  is  called  whenever  the  calling  program  requires 
the  fractions  burned,  the  Z  array,  of  the  M  propellants  being  burned.  A 
glance  at  the  listing  reveals  that  FORMT  is  essentially  one  large  loop, 
processing  each  of  the  M  propellants  in  turn.  All  computations  for  the 
Jth  propellant  arc  bypassed  if  either  it  is  entirely  consumed,  that  is, 

Z(J)  i  1 

or  if  logic  elsewhere  in  IBHVG  dictates  that  the  condition  flag 

ICFLAG(J)  -  0 

which  signifies  that  the  ignition  criterion  for  the  Jth  propellant  has 
not  been  met.  However,  if  the  propellant  has  ignited  but  is  not  yet 
completely  burned,  execution  continues.  The  depth  burned  for  propellant 
J  is  extracted  from  the  Y  array,  having  been  calculated  elsewhere  prior 
to  entrance  to  FORMT,  and  twice  this  value  is  stored  in  U.  Next,  N(J), 
the  grain  code  for  propellant  J  in  Table  1  is  examined,  and  a  branch 
taken  to  the  appropriate  area  of  FORMT,  There,  the  fraction  burned, 

Z(J) ,  is  determined  based  on  U,  initial  grain  dimensions,  and  various 
control  variables,  followed  by  a  jump  to  the  bottom,  of  the  loop.  After 
all  M  propellants  have  been  processed,  a  return  is  made  from  FORMT  to 
the  calling  routine. 

As  Table  1  notes,  there  are  ten  different  grain  types  dealt  with  1 
by  IBHVG.  The  remainder  of  the  report  will  be  oncemed  with  consider¬ 
ing  the  types  in  turn,  devising  a  method  of  £i  ling  Z  for  them  at  any 
instant,  and  developing  the  FORMT  coding  based  on  these  calculations. 


III.  GRAIN  TYFB1:  SBVBN-PERFORATED  CYLINDER 

An  end  view  of  this  grain  appears  in  Figure  1.  It  is  a  cylinder 
pierced  longitudinally  by  seven  cylindrical  holes  arranged  at  the  center 


TABLE  1.  IBHVG  GRAIN  GEOMETRIES  AND  CODES 


GRAIN  CODE 


GEOMETRY 

7 -perforated  cylinder 
1 -perforated  cylinder 
cord 

rectangular  strip 
sphere 

slotted  tube 
37-perf orated  hexagonal 
19-perf orated  hexagonal 
19-perforated  cylinder 
7 -perforated  hexagonal 


and  vertices  of  a  regular  hexagon.  Initially,  suppose  the  grain  to  have 
length  GL(J) ,  diameter  D(J),  perforation  diameter  PD(J),  inner  web 
WI(J),  and  outer  web  WO(J) . 

On  the  first  call  to  FORMT,  i.e.,  at  the  beginning  of  the  first 
time  step,  COW-ON  variable  1  is  set  to  0.  On  succeeding  passes,  I  will 
be  greater  than  zero,  thereby  forcing  a  branch  to  line  SO,  bypassing  the 
initialization  section  which  is  described  next. 

First,  the  middle  web,  WM(J)  is  zeroed  as  it  is  not  used  by  this 
grain  type.  Then  the  diameter  is  recomputed  to  absorb  any  inconsistencies 
in  grain  dimensions.  Figure  1  makes  the  calculation 


D(Jf)  «  3  •  PD(J)  +  2  •  (WI(J)  +  W0(J)), 


clear.  Next,  the  initial  grain  volume  is  determined  and  stored  in 
VO(J) .  The  end  area,  B,  of  the  grain  is  given  by 


B  «  ir/4  •  D(J)2  -  7  •  it/4  •  PD(J)2, 


so  the  initial  volume  is 


V0(J)  -  GL(J)  •  E 


-  .25  *  it  «  GL(J)  •  (D(J)2  -  7  •  PD(J)2). 


Note  that  we  assume  ir  is  passed  through  COMMON  Variable  PI. 

The  method  developed  in  Reference  4  is  employed  and,  in  what  follows, 
its  notation  and  nomenclature  are  used  freely.  For  purposes  of  calcu¬ 
lating  the  end-area  of  the  grain,  it  is  partitioned  into  one  class  of 


'1  "T» y  1  ««.*.,  -n-r-...-, 


outer  slivers  congruent  to  ABDC  and  one  class  of  inner  slivers  congruent 
to  OAB. 

Array  I  FLAG  is  the  branch-flag  array  and  two  positions  are  ini¬ 
tialized  to  -1  for  this*  the  Jth  propellant. 

Array  SF  is  the  "side"  array.  Letting 


WW  »  WI(J)  +  PD(J) 


(4) 


represent  the  distance  between  centers  of  adjacent  perforations*  all 
three  sides  of  each  sliver  in  both  classes  are  set  to  length  NW. 

The  last  calculation  in  the  initialization  section  is  the  deter¬ 
mination  of  the  grain’s  web,  the  shortest  distance  through  the  grain  to 
be  broached  as  the  grain  bums.  This  is  given  by 


WBBC(J)  -  min  (WO(J) ,  WI(J),  GL(<J))* 


(55 


the  length  being  included  to  account  for  very  short  grains. 

The  post-initialization  coding  starts  at  line  50.  The  reduced 
grain  length*  reduced  grain  diameter*  and  increased  perforation  diameter 
are 


GRL  -  max  (GL(J5  -  U,  05 

OD  -  D-U  (6) 

PRFD  -  PO+U, 


respectively,  the  maximum  being  taken  to  guarantee  a  non-negative  grain 
length.  If 


0  <  U  <  WEBC(J), 


the  web  is  not  broached  and  we  determine  ond-area  of  the  grain*  B,  by 
the  same  method  as  in  the  initialization  section,  so 

E  -  ir/4  •  OD2  -  7  •  r/4  •  PRFD2  (7) 

yielding 

Z(J)  -  1  -  GRL  •  E/VOCJ).  (8) 
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However,  if 


WEBC(J)  <  U, 


we  branch  to  line  60  and  employ  routines  GENIS  and  GENOS  as  described  in 
Reference  4  with  one  minor  modification.  Since  IBHVG  does  not  use  grain 
surface  area,  the  argument  to  return  that  quantity  is  not  present  in  the 
call  and  the  coding  in  GENIS  and  GENOS  which  yields  sliver  surface  area 
is  commented  out  as  can  be  seen  from  APPENDIX  A,  The  AW  array  is  the 
work  array  and  sliver  volumes  for  the  inner  and  outer  slivers  of  the  Jth 
propellant  are  returned  through  the  appropriate  cells  in  array  GV. 

Total  sliver  volume  of  a  given  class  is  six  times  the  volume  of  each 
sliver,  yielding 


Z(J)  -  1  -  6  •  (GV(1,J)  +  GV(2,J))/VO(J). 


O) 


A  branch  is  then  made  to  the  bottom  of  the  loop. 


IV.  GRAIN-TYPE  2;  SINGLE- PERFORATED  CYLINDER 

Initially,  this  grain  is  a  cylinder  of  length  GL(J)  and  diameter 
D(J)  pierced  by  a  cylindrical  hole  of  diameter  PD(J) ,  the  centers  of  the 
two  cylinders  coinciding.  Let  the  web  be  WI(J). 

As  with  the  case  of  the  seven-perforated  cylinder,  an  initialization 
section  is  processed  if  I  «  0.  The  middle  web,  WM(J) ,  and  outer  web, 
W0(J) ,  are  zeroed  as  they  are  unused  and  the  diameter  recalculated  as 

D(J)  -  PD(J)  +  2  •  Wt(J),  (10) 


to  absorb  any  error  in  grain  dimensions.  Then 


WEBC(J)  *  min  (W1 (J) ,  GL(J)) 


(11) 


is  twice  the  depth  burned  at  grain  extinction.  Letting  E  represent  the 
initial  end  area  of  the  grain,  we  have 


B  ■  ir/4  •  D(J)2  -  tt/4  •  PD(J)2, 


(12) 


and  the  initial  volume  is 


VO(J)  -  GL(J)  •  E 


-  fl/4  •  Gl(J)  •  (D(J)2  -  PD(J)2). 


(13) 


The  factor  it/ 4  is  not  present  in  the  coding  as  it  will  cancel  in  the 
calculation  of  Z(J). 

The  post- initialization  coding  starts  at  line  190.  The  end-area, 
B,  of  the  grain  burned  to  a  depth  of  1/2  U  must  be  determined.  First, 
it  is  cleared  by 


E  -  0. 


If 

0  <  U  <  WBBC{J) , 


then  we  override  this  value  of  zero  for  E.  The  reduced  grain  length, 
reduced  grain  diameter,  and  increased  perforation  diameter  are 


GL(J)-U 

D(J)-U 

PD(J)+U, 


respectively,  so,  as  above, 


E  -  ir/4  .  (D(J)  -  U)2  -  ff/4  •  (PD(J)  +  U)2 
-  it/ 4  •  ( (D(J)  -  U)2  -  (PD(J  )  ♦  U)2), 


and  fraction  burned  is 


(14) 


Z(J)  -  1  -  B  •  (GL (J)  -  U)/V0(J). 


(15) 


Note  that  the  factor  n/4  in  the  calculation  of  B  is  absent  in  the  coding 
as  it  would  cancel  in  the  Z(J)  calculation. 


V.  GRAIN  TYPE  3:  CORD 

Initially,  this  grain  is  a  solid  cylinder  with  length  GL(J)  and 
diameter  D(J) . 
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If  COMMON  variable  I“0,  indicating  the  first  call  to  FORMT,  an 
initialization  section  is  processed.  The  perforation  diameter  and 
inner,  middle,  and  outer  webs  are  zeroed  for  this,  the  Jth  propellant, 
as  these  quantities  are  unused.  The  initial  volume 

V0(J)  -  GL(J)  •  ir/4  •  D(J)?  (16) 


is  that  of  a  right  circular  cylinder.  In  the  coding,  the  tt/4  is  omitted 
because  it  will  cancel  in  the  calculation  of  Z(J). 

Post-initialization  code  starts  at  line  220.  The  reduced  grain 
dimensions 


GL  (J)-l' 
D(J)-U 


are  used  to  produce  the  fraction  burned 

Z(J)  «  1  -  max(GL(J)-U,0)  •  max(D(J)-U,0) 2/V0(J)  (17) 


where  th«i  tt/4  factor  is  absent,  having  cancelled  with  a  like  factor  in 
V0(J). 


VI.  GRAIN  TYPB  4:  RECTANGULAR  STRIP 

Initially,  this  grain  is  a  solid  rectangular  parallelopiped  with 
length  GL(J)  width  I)(J),  and  depth  WI(J). 

As  above,  variable  1=0  signals  processing  of  an  initialization 
section.  First,  all  grain  dimensions  not  used  are  zeroed  and  the  initial 
grain  volume 


VO(J)  *  GL(J)  •  D(J)  •  WI(J),  (18) 


is  determined. 

The  post-initialization  code  starts  at  line  440.  The  reduced  grain 
dimensions,  namely. 


GL(J) -U 
D(J)-U 
WI(J)-U, 
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are  multiplied  together  to  produce  the  volume.  And  thus  the  fraction 
burned  is 


Z(J)  =  1  -  max(GL(J)-U)  •  max(D(J)-U)  •  max(WI(J) -U)/VO(J)  (19) 

where  the  maxima  are  taken  with  0  to  guarantee  a  non-negative  grain 
volume . 

VII.  GRAIN  TYPE  St  SPHERE 

Initially,  this  grain  is  a  solid  sphere  of  diameter  D(J) . 

On  tho  initialization  pass,  signaled  by  I  «  0,  all  grain  dimensions 
except  grain  diameter  are  zeroed  because  they  are  unused. 

The  reduced  grain  diameter  is 
D(J) -U. 

Thus  fraction  burned  is 

Z(J)  -  1  -  1/6  •  ir  •  max(D(J)  -  U,0)X/(l/6  •  ir  •  D(J)3)  (20) 

-  1  -  (max(D(J)  -  U,0)/D(J))3, 

where  the  maximum  is  taken  to  guarantee  a  non-negative  reduced  grain 
diameter. 


VIII.  GRAIN  TYPE  6:  SLOTTED  TUBE 


An  end  view  of  this  grain  appears  in  Figure  2.  It  is  an  ordinary 
single-perforated  grain  with  a  slot,  the  sides  of  which  are  parallel. 
Initially,  suppose  the  grain  to  have  length  GL(J),  diameter  D(J) ,  web 
WI(J),  perforation  diameter  PD(J) ,  and  slot  width  WM(J) . 

As  with  the  other  grain  types,  an  initialization  section  is  processed 
if  COMMON  variable  I  »  0.  First,  the  unused  dimension  W0(J)  zeroed. 

Let 


PRFD  -  .5  •  WM(J) ,  (21) 


be  half  the  slot  width. 


PRFD 


SO  »  .5  •  D(J) 


(22) 

be  the  grain  radius,  and 

SI  ■  .5  •  PD(J) ,  (23) 

be  the  perforation  radius .  In  right  triangles  OBD  and  OAC  we  have 
9  -  tr/2  -  acos  (PRPD/SO) 

(24) 

a  ■  ir/2  -  acos(PRFD/SI) , 

respectively.  As  before,  the  code  as  suae  s  PI  contains  ir.  Denote  by  B 
the  end-area  of  the  grain;  then 

B  *  area  reflex  sector  BOP 

-  area  reflex  sector  AOE 

-  2  *  area  triangle  AOB 

-  1/2  •  (2ir  -  28)  •  SO2  -  1/2  ♦  (2ir  -  2a)  •  SI2 

-  2  •  1/2  •  PRFD  •  (SO  •  CQS(6)  -  SI  •  cos(a)) 

*  (ir  -  6)  •  SO2  -  (it  -  a)  •  SI2 

-  PRFD  •  (SO  •  cos (8)  -  SI  •  cos(a)) .  (25) 

Thus  initial  volume  is 

VO(J)  »  GL(J)  •  B.  (26) 

Finally, 

WBBC(J)  ■  ain(GL(J) ,  WI(J)),  (27) 

is  twice  the  depth  burned  at  grain  extinction. 
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The  post-initialization  section  of  coding  starts  at  line  S10.  The 
end-area  of  the  grain,  E,  is  cleared  by 

B  ■  0. 

If 

WBBC(J)  <  U, 

then  the  grain  is  extinguished  and  a  branch  is  taken  to  line  520, 
producing  a  fraction  burned  of  1.  Otherwise,  let 

PRFD  ■  .5  *  (WM(J)  +  U),  (28) 

be  half  the  new  slot  width, 

50  -  .5  •  (D(J)  -  U),  (29) 

be  reduced  grain  radius, 

51  -  .5  •  (PD(J)  +  U)  (30) 

be  perforation  radius,  and 

GRL  -  GL(J)  -  U,  (31) 

be  reduced  grain  length.  Then  6,  a,  and  E  are  calculated  as  done  in  the 
initialization  section  and  fraction  burned  is 

Z(J)  -  1  -  B  ♦  GRL/ VO (J) ,  (32) 

IX.  GRAIN  TYPES  7,  8,  AND  10:  R0UNDBD-HBXAGONAL 

Figure  3  depicts  one-sixth  of  a  rounded-hexagonal  grain  with 
nineteen  perforations.  Let  grain  length  be  GL(J),  all  inner  webs  be 
WI(J),  the  outer  web  W0(J),  and  perforation  diaaeter  PD(J).  Note  that 
the  rounded  "comers"  are  bounded  by  inside  and  outside  circular  arcs  of 
radii  1/2  *  PD(J)  and  WO(J)  +  1/2  *  PD(J),  respectively,  as  shown  in 
sector  CDB. 


Figure  3.  One-Sixth  of  Hex  Grain  -  End  View 
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The  7-,  19-  and  37-perforsted  hexagonal  grains  are  but  special 
■embers  of  a  family.  By  induction  on  K,  the  number  of  hexagonal  rings 
of  perforations,  we  have 

NP  ■  number  of  perforations  ■  3K2„+  3K  +  1  (33) 

SI  >  number  of  inner  slivers  ■  6a 
SO  ■  number  of  outer  sliver#  ■  6K 
nusber  of  comer  slivers  »  12, 


for  K  ^  0.  The  inner  slivers  are  all  congruent  to  the  area  within  ABC, 
the  outer  slivers  to  the  area  within  BCBF,  and  the  comers  to  CDB. 

Routine  FORMT  currently  handles  cases  K  ■  1,  2,  and  3  yielding,  respectively, 
NP  •  7,  19,  and  37,  SI  -  6,  24,  and  54,  and. SO  -  6,  12,  and  18. 

Entry  for  hex  grains  is  made  at  statements  240,  250,  or  600  in 
FORMT  where  SO  and  SI  are  set  to  the  appropriate  values  and  a  branch 
taken  to  line  260.  Here  WW  is  the  distance  between  centers  of  adjacent 
perforations  and  NW2  is  its  square,  PRFD  is  the  increased  perforation 
diameter  and  PRPD2  its  square,  GRL  is  the  reduced  grain  length,  snd  NOD 
is  the  outer  web.  The  end  area  of  the  grain,  B,  is  initially  cleared  to 
zero. 

Figure  4  shows  an  inner  sliver  during  burning.  By  definition  of 
cosines, 


6  «  2  •  acos(min(l,  CT/C7))  (34) 

-  2  *  acos(min(l,  WW/PRFD)), 


where  the  minimum  is  taken  with  1  in  case  0  is  degenerate.  Figure  5 
shows  an  outer  sliver  burning.  By  definition  of  cosines, 


a  ■  acos(min(l,ECy£H))  (35) 

-  acos(min(l,(PD(J)  +  2  •  N0D-U)/PRFD|) , 


where  the  minimum  is  taken  with  1  in  case  ot  is  degenerate.  These  angles 
are  sufficient  to  calculate  E. 

Consider  the  twelve  comers.  If 
NOD  i  U, 


they  are  completely  consumed  and  their  contribution  to  E  is  zero. 
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Figure  4.  Hi 


-i:,. 


Otherwise,  by  subtracting  the  area  of  concentric  circular  sectors,  and 
recalling  that  it  is  in  PI, 

E  -  12  *  (1/2  •  u/6  •  (1/2  •  PD(J)  ♦  WOD  -  1/2  •  U)2  (36) 

-  1/2  •  tr/6  •  (PRFD/2)2 

■  .25  •  it  •  ( (PD(J)  +  2  •  WOD  -  U)2  -  PRPD2). 

Consider  now  the  inner  slivers.  A  DATA  statement  provides  values 
of  square  root  of  3  in  RT  and  it/ 3  in  PI3,  variables  required  for  the 
remainder  of  this  analysis.  If 

ir/3  <_  e, 

then  all  inner  slivers  are  consumed  so  a  branch  is  taken  to  line  270. 
Otherwise,  the  area  contribution  of  the  SI  inner  slivers  must  be  added 
to  B.  Hence,  referring  to  Figure  4, 

B  ■  B  +  SI  •  (area  inner  sliver) 

■  E  +  SI  •  (area  traingle  ABC  -  3  •  area  triangle  JBC 

-  3  ♦  area  sector  JCK) 

■  B  +  SI  *  (1/2  •  WW  •  (1/2  WW  ^ 

-  3  •  1/2  •  WW  •  (1/2  •  PRFD  •  sin  (6/2)) 

-  3  •  1/2  •  (tt/3  -  6)  *  (PRFD/2)2). 

Then,  by  substitution  in  the  line  with  sin  (6/2)  using 

WW  -  PRFD  •  cos (6/ 2) ,  (38) 

and  recalling  the  half-angle  formula  for  sin, 

B  -  B  +  SI  •  .25  •  (WW2  1.5  •  PRFD2  •  (sin(6)  +  ir/3  -  6)).  (39) 


Lastly,  consider  the  outer  slivers.  When 


or,  equivalently 


a  1/2  o-e) , 

.  I 


then  the  outer  slivers  are  extinguished  and  a  branch  is  taken  to  line 
280.  Otherwise,  the  area  contribution  of  the  SO  o  ‘*r  slivers  oust  be 
added  to  E.  Hence, 


E  ■  B  +  SO  •  (area  outer  sliver)  (40) 

»  B  t  SO  •  (area  rectangle  BCBF  -  area  triangle  BCP 

-  2  •  area  triangle  MEC  -  2  •  area  sector  MCP) 

-  E  +  SO  •  (WW  •  (PD(J)  +  2  •  W0D-U)/2 

-  1/2  .  WW  •  (1/2  •  PRFD  •  sin  (0/2)) 

-  2  •  1/2  •  (PD(J)  +  2  •  WOD-U)/2  .  (1/2  •  PRFD  •  sin(o)) 

-  2  •  1/2  •  (it/ 2  -  a  -  9/2)  •  (PRPD/2)2) 

-  B  +  SO  .  .125  •  (2  •  (2  •  WOD  ♦  PD(J)-U)  •  (2  •  WW  -  PRFD  •  sin(a)) 

-  PRFD2  •  (sin(0)  ♦  *  -  2  •  a  -  B))« 

This  concludes  the  calculation  of  E.  If  variable  I  >  0,  then  the 
initialization  section  is  bypassed  and  Z  *  1  -  E  •  GRL/V0(J).  Otherwise, 
the  grain  diameter  D(J)  is  recalculated  for  printout  purposes,  As  in 
the  case  of  the  other  grain  types,  the  initial  grain  volume  must  be 
computed , 

Initially,  we  have 


a*0 

9-0 

PRFD  ■  PD(J) , 


and  all  slivers  and  corners  are  present.  This  implies 

E  -  .25  •  tt  •  ((  PD(J)  +  2  •  WOD)2  -  PD(J)2)  (41) 

+  SI  •  .25  •  (WW2  •V7-  1.5  •  PD(J)2  •  t/S) 

+  SO  -  .125  •  (2  •  (2  •  WOD  +  PD(J))  •  2  •  WW  -  PD(,J)2  •  w) 


-  w  •  WOD  •  (WOD  +  PD(J) )  +  .25  •  SI  •  WW2  •  V? 

-  .125  •  *  •  (SO  +  SI)  •  PD(J)2  +  SO  •  WW  •  (WOD  +  1/2  PD(J)>, 


and 

VO(J)  -  GL(J)  •  E.  (42) 


Finally,  a  check  mist  be  aade  to  see  whether  the  outer  sliver 
remains  within  its  original  bounds,  for  if  it  does  not,  the  above  analysis 
is  invalid.  That  is,  the  point  P  in  Figure  5  mist  not  reach  the  midpoint 
of  6?  before  tho  perforation  centered  at  C  does.  Equivalently,  we  mist 
have 


!(T  >  1/2  CF  -  1/2  PD(J)  (45) 

WOD  +  1/2  PD(J)  >  1/2  WI  (J) 

2  *  WOD  +  PD(J)  >  WI(J), 


and  conversely.  Otherwise,  FOMfT  will  stop  with  an  error  message. 


X.  GRAIN  TYPE  9}  NINBTBBN-PBRPORATBD  CYLINDER 

Initially,  this  grain  is  a  cylinder  pierced  by  nineteen  uniform 
cylindrical  perforations.  An  end  view  is  shown  in  Figure  6.  Its  length 
is  GL(J) ,  diameter  D(J),  and  perforation  diameter  PD(J).  It  is  charac¬ 
terized  by  three  webs:  inner,  WI(J);  middle,  WM(J) ;  and  outer,  WO(J) . 

If  COMMON  variable  I  is  greater  than  zero,  the  initialization 
section  is  bypassed  by  a  branch  to  line  320.  Otherwise,  the  grain 
diameter  is  recalculated  to  absorb  any  inconsistencies  as 


D(J)  «  S  •  PD(J)  +  2  •  (WI(J)  +  WM(J)  ♦  W0(J)),  (44) 

and  the  grain  volume  is  computed  by 

V0(J)  ■  GL(J)  •  (ir/4  •  D(J)2  -  19  •  ir/4  •  PD(J)2)  (45) 

-  .25  •  it  *  GL(J)  •  (D(J)2  -  19  •  PD(J)2). 
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The  r attainder  of  the  initialization  section  and  the  whole  of  the 
post-initialization  coding  is  drawn  directly  from  the  worked  example. 
Section  IV,  in  Reference  4  with  IFLAG,  SF,  AW,  GV  being  the  flag,  side, 
work,  and  volume  arrays.  The  only  novelties  are  the  addition  of  an 
additional  subscript  for  the  Jth  propellant  and  the  removal  of  the  sur¬ 
face  area  computation.  The  latter  was  discussed  in  the  case  of  the 
seven-perforated  grain,  above. 
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APPENDIX  A 


LISTING  OF  FORM-FUNCTION  SUBROUTINES 


FOR  Mr  <007  Z9-APR-8Z 

SUBROUTINE  FOUNT 
COINON/IBNITE/ICFLADMI 

CONNDN/RUNBLK/XCi ZQt* PR(ZO)*H(9)*8TV(9)«CO)* 
9Fm.TO(9»*ALP(9>*BET(9)*8At9>*COV«9>*RMO(9>» 
teL<9).Dm*P0<9)*tfZ<9)*NNt9>*U0t  9 1. 81*62(21* 

Mitt#  90*  T  W*  D8*DL*GLR*  CV*DT*  MP*NFR*N*L V3PK* ICQOE t9 1* 
ITHRESH(9)*IIRMO(9l*BRHP(9*9)*BRHR(9*9)tAVBWEB(9) 
C0NH9N/XTRA/PT(2>»PB(2t,XLST*PRLST,DEN2*N(2)»RES.A8N» 
lt»CT* AF*PN*P6R(2)*0EN1*HCL* AC  ( 9  )*  BC  (5 ) *CC  1 9  )* DC  (9 1* 

BIC(  9)  »  ZI9 1*  P0V*FH1*  PT  PF*  P6PP*  HPF*  PBRPP#TPF»  VPP*  XPP*  AFP* 
BAPK*Z1PP*Z2PP*Z3PP»ZAPP*Z9PP*  COEFF*DEN3*CTN*PI*VO(9l* 
IXR.NN.VtBl.DYtBt 

DIMENS  ION  WEBCOl.SFt  9.  A*  9 1  *  AWU*  A*  5)  ,  IFL  A8(A,  91. 69(  4*  91 
DATA  RT 71*  7 B20908 OB/*  P 13/ 1*04  7197991/ 

FIND  FRACTIONS  BURNED  IU)  FOR  1*1*2*  ...*N  PROPELLANTS 
00  1900  J  *  1  *N 

B7PASS  IF  JTH  PROPELLANT  UNI8N1TED  OR  EXTIN6UISHED 
IFtZt  J».&E.1..0R.ICFU8U).E1.0)  SOTO  1000 

SET  U  ■  2  X  DEPTH  BURNED  (  BRANCH  ON  8RAIN  TYPE 
U*2.»7(JF3) 

SOTO  (20*190*  210*  ISO*  490*  900*  290* ZAO*  110*6001  *N(J  > 

CODE  li  7-PERF  8RAIN 
IF ( I, ST .0  )  OOTO  90 
NN(  Jl-O. 

OUt»9,9PO(  J)A2.*(MI(1)»N0(  J)  I 

Y0< J»..29*PI*OL (119(0(1 »**2-7. APDI1)942) 

tFLA8ll*Jt--l 

lFLAJt2.JI*-l 

NU-WIU)FPD(J) 

DO  JO  K»l*  9 
$F(K*1*J)*MN 
SF(K.2*1I*WW 
CONTINUE 

1EBC(J)*AMINL(M0(1)*XI(1)*QL(JH 

9RL«ANAU(SL(J)-U*0.I 

3  0  ■  0  (  4 1  —II 

FRFO-POtmU 

|F(U.OT*tfEBC(Jt)  SOTO  60 

E •«  2JAPI* (00**2-7*9PRF09*2l 

Z( JI«1.-6RL*E/VU(1> 

OOTO  1000 

CALL  OFNlS(SFU,l,J»,A'm*l*J>*PRFO*ORl#IFLAO(l»J>*8V(l*lH 
CALL  OEMOKSFU*  2*1).  AW(1*2»J  I.PPFD.8RL* .9*00* IFLA8 < 2* J I »0V( 2.1 ) ) 
l(  il*l.-6.«(8V(l*  JIA8VI2*  J)  1  /  VO  ( J  I 
OOTO  1000 

CODE  Zl  l-PERF  DRAIN 
IFU.ST.o)  SOTO  190 
tMdl'O. 

V0(  41  *U* 

0m.p0tm2.9wn  n 

WE8C(  1>*AH1NI(GL(  Jl.WHJH 


’vTr?7jP'!r~:  p  «c  ' , 


no 


c 

CM* 

no 


IZO 

C 

CM* 

2*0 


C 

CM* 

290 

260 


no 


290 


290 

c 

CM* 

Xlo 


219 


90(41  *61(4  l*(0(4)M2-ROtilM2  I 
EM. 

lMU.tE.N«*C(4))  l*tO(JI-U)Mt-f*Ot41*U)**« 

zomi.-emsiui-umnoui 
SOTO  1000 


caoe  i<  curd  saaih 

IXU.ST.O)  SOTO  220 
*0(41*0. 

91(41*0. 

N0(J)*0. 

HNt J)*0. 

VOt  41  *91(41*0(4 1**2 

1 1 4  1*1  .-ANAX1  IBl(  41-4*9*  )♦( ANAX1 (Ot4 1-0*0. 1  >•♦* # VO (49 
SOTO  1000 

cool  S'  ROUND-HEX  19-FEM  0*1  IN 

so-n. 

*1*16. 

SOTO  260 

COOE  7 1  ROUND-HEX  17-RIRf  SKAIN 
10*19. 

St*94» 

NH*NI  (41**0(41 

nm*nn**2 

RR*0*RD(4)»U 

ARF02MM0M2 

SRl*ANAXl(Sl(4l-0*0.l 

400*110(4) 

£•0. 

rM«U-t.MC0S(ANtNl(4W/**90,l.)) 

4L  RNA* AC99 1  AHtNl ( (2.  *MOO*  *D  (4)-U)/A**D#l.l) 
l F  (U« IT. 400 )  6*.2J**I*(<2.«NJO**Ot4l-U)M2“MFOtl 
IF<MFTA.6E.*tll  SOTO  270 

«.!*JI*.25*(WW2*RT-1.  J**RFD2*<SlM(74iTA)**l»-TM«TAn 
I r ( At. RHA.SE •  .9*(  PI-THETA))  SOTO  111 

E*E*J  9*«129*t  2.*(  2.*NOO**D(  4  )'*U)*tI  •♦NN-**RD*SIH(  ALMAl  1-MF02* 
ItStNITHStAI  ♦*I-2.*Al*HA-?H2T* ) ) 
tAU.8T.9l  SOTO  290 
NN( 4) • MI  (  4 ) 

tRIHI41.E0.10)  0(4)«).**OI4)*2.*H1I4)*2.*MOD 
IR(N(4I.E9.S)  0( 4 )*9. 4*0(41 *4 ,*NI (4 ) *2.*W00 
IF(HU).EO.T)  OI4)*7.MOt4)*6.4Ml(4t*2.*MOD 
V  0 ( 4 1 *E*ftL ( 4 ) 

I F ( 2. *WOD 4*0(4). OE. Nl (4 ) )  SOTO  190 
CALL  8A«HSS(**F0RHT4  BAO  HEX  MQA'I 
2(41*1. -f  *SKIM14> 
son  looo 


COSE  91  19-RERF  SRAH 

IFU.Ol.O)  SOTO  120 

no  l-9.»PO(41*2.«  (Nil  4 )*MHl  4 )  **0(4)  > 

NOt  J)..*5**14SU4)*(0(4l**2-19.**0t4)4*2) 
00  119  K»  l»  4 
t *L AS (K« 4 )«-l 

runm 
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XMl.lt4)*XI(J)*fD(4> 

tMMiJI-lfllil.J) 

JMIr  1»4)*SM1.1.4I 

SMt»t*4)*SMl»«»  J) 

SM)»I»4I*SM1»1*4I 

SMliB.4l*RD(4>«.94<«t(4)4iH(4)) 

JMt.  I»4I-5M1»I»4> 

SMI.  I»4)*WN(4)4RD<4) 

SM1.4.4I*SM1>S>4) 

SMt.4.4l*1.4$Hli>»4) 

4M).4.4)*SMl.t.4)4RT 

xiicuj.AMNMXouuxiiui.Mut.jMi.iui-HUJ.iMi.tut-pouj. 

muii 

Ito  BAL*ANAXl(0l(41-l)p0.) 
oo*o(4>-u 
MRD*R0(4)4U 

IMU.8T.XIICU))  «0TU  ISO 

l<4)*&.-0RL4R/V0(4l 
SOTO  XOOO 

no  on  no  k« i#  i 

CALI  OfiMIS<SMlfK.4)»AWtl.X»4>.MI:O.ORLimAO(K.4>.BV(K»4l> 

140  CONUHUI 

CALL  atNOI(SMl#4.4).  Am»4.4)»MPD»m».J*00,IRLABI4.  JltlTUtDI 
l(J>.l.-(b.*(6V(1.4l  +  6V(t«J)) *11 <4(01(1.4)401 (4.4)1 1/10(4) 

6010  XOOO 

c 

C444  coot  41  RICTANOULAR  STRIP  IRA  IN 
410  1MI.0T.0)  10 TO  440 

90(41*0. 

10(41*0. 

11(41*0. 

101 4) *6L( 4  >*0(4)411(4  I 

440  im«l.-AHAXl(BL(4)-i4#0.>4ANAXl(0(4)-U*0.  )* 

*44*Xi(1t(4)-U»0.)/10(4> 

QOTO  1000 

c 

:♦♦♦  CODI  91  SMCNICAL  BRAIN 
410  IMl.OT.OI  OOTO  440 
BLUI  *0. 

9  01 4)  *0. 

11(4)*0. 

40(4)  *0. 

44(41*0. 

440  IU>*l.-(4nAXi<0<4)-U.0. 1/0(4)  >44} 

OOTO  1000 

m 

C*4*  C 301  4)  ILOTUO-TOBI  BRAIN 
90)  (Ml.  11.0)  4010  9X0 
00(41*0. 

PA90*»»*41I4I 

10«.9*0(4> 

91*. 1490(4) 

H«TA«.5<‘Pl-*COS(PAItO/SO> 

AL91A.,»*Pl-AeOS(9R9D/Sl> 

10(41*61(4 I ♦( (9 1»1HIT  A)4S04*I-(9I”ALPHAI*$I *4t» 


»MrD*<»o*comn«ni-si*cjiuiMAm 

tfiicm'AnmiuutHiun 

no  i*o. 

!Mu.«T<4iicun  ton  no 
MP0*»>*(NN(JI4U) 

(0««S4(9( JI"W> 

n«.f*(»om«ui 

ONlMim-U 

WITA'tlAf  l-ACDSIMNO/IQI 

I-AlMAimMl-fMOMSO* 
ICOI (THtTM-SIKOt UlMUH 
*19  im*l»-«*»U/V0(9) 

•ora  iooo 

c 

CAM  COOI  101  OOUMO-MIK  7-OIIP  HUN 
400  $0*4. 

Sl*4. 

«oro  too 

c 

1001  G04T1HUC 
NITUNN 
(NO 


3S 


rto n  non  on n  no o  ono  oo n  n  ooooo  n  nnnnn 


C  tl1IS.»01  11-0IC-7S 

SUSR9UT1NS  SCNIS(S»A*PRPI>*IRL*  IPlAS.T'OL) 

SUSROUTINI  OSINISOi  CALCULATI  SURPACI  (IIA  AN0  VOlUNt  PM  A 
•INIRAL  INMI  SLIVIR  OP  A  IUM1M  MAIN 
NIYN  LINQTri  •  ML  A  PIPP  01AN  •  PSP*. 

DININSION  SII)*AtA) 

DATA  PUP  l*IT07f Al  l 

IPttPLAA)  10.70*10 

INITIAL  PASH  IPIA0  MAS  SI  I  NIIATIVI  0V  CALL  INI  ROUTINI. 

STOPS  ANSLIS  A1.A7.AI  AN0  ARIA  OP  TR1AMLI 
WITH  SIDIS  I  (Hi  St  t)*S(S)  INTO  AI1)m..iA((I 

10  Ain  •  icQSttst<)«*t+s(s>**t-s<iiMimt.*stti*s(>m 
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•  ••CHICK  POR  I  PROP  CONDITION!  PINO  IP  TRIAN1LI  ACCIPTASLI... 
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00  10  1  ■  1*S 
IPIAID.LT*  .SAPID  J  ■  JPl 
10  C0HTINUI 

tPli.QT.il  CALL  IRRHSSDPQINIS*  IRR0R  «) 

tP  OK*  SIT  PLAO  10  IfPO  TO  SVPAIS  INITIALIIATION  rilRlAPTIR. 

I  PL  AS  •  0 

S ICCEIDINQ  PASSU  UNTIL  lURNOUTl  PINO  AUXILIARY  ANSLIS 

10  TAU12  -  ACOSI ANIN1I1*  *3  It )VPRP0) I 
TAU11  ■  ACOSI AHINlll. *  Sil >/PRP0) ) 

TAU7S  •  ACOSIANINlll.iSID/PRPOH 

...AND  OKANCH  IP  SlIVIR  PAILS  SURNOUT  CRITBRIA. 

tPITAU12PTAUllPTAUIl.LT.PI2  .AND*  ORL.QT.O.)  SOTO  IS 
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tPLAS  •  1 
SOTO  10 
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10  I  ■  Al A)“*IS*PRPO*(SI 1)PSINITAUIS)PSI2 I*S1NITAU1I) 
tPl(l) PS INI T  AU1I »P PR POP  I  PI P"TAU1I-TAUIS~TAUIS) I 

VOL  •  IPBRL 
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C  SUSP  •  •••IpQRLPPRPO* IPII-TaUII-TAUIS-TAUII ) 
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subroutine  hnos(s.a.errd.srl*rad>ieiai.vol) 

SUBROUTINE  «SIN0**l  CALCULATE  SURFACE  ARIA  AND  VOLUME  E0R  A 
tEHIItl  OUtIR  SLIVER  QE  A  BURNIMS  MAIN 
WITH  LINRT H  -IRLf  RADIUS  •  RAO.  AMO 
EIRE  01 AM  •  EREO. 

DIMERS  ION  Stt)»At4> 

IE  I  IE  LAI)  10.M.S0 

INITIAL  EASSI  IELA0  HAS  SIT  Ml  DATIVE  IT  CALL  IMS  ROUTINE. 

STORI  ANBLIS  A1.AE.AI  AND  ARIA  OR  TRIAMf LI 
HI TM  SIDES  SUI.Sttl.Stl>  INTO  At  LI..... A(4I 
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1ELAS  -  1 
tUTO  SO 
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LIST  OF  SYMBOLS 


Symbol 

Array 

Description 

M 

N  (NO) 

Number  of  propellants 

Z 

Y  (YBS) 

Fractions  burned 

ICFLAG 

Y 

Ignition  threshold  flags 

N 

Y 

Grain  codes 

U 

N 

Twice  depth  burned 

I 

N 

Flag,  equals  0  on  initial  pass 

VO 

Y 

Initial  grain  volumes 

IFLAG 

Y 

Slivering  flags 

GV 

Y 

Sliver  volumes 

WBBC 

Y 

Overall  grain  webs 

SF 

Y! 

Work  arrays 

AN 

y) 

RT 

N 

3 

PI 

N 

PIS 

N 

ir/3 

GL 

Y 

Grain  lengths:  all  codes  except  5 

D 

Y 

j  Grain  diameters:  all  codes  except  4 
l  Grain  widths:  code  4 

PD 

Y 

Perforation  diameters:  codes  1,2, 6, 7, 1,9 

NI 

Y 

j  Inner  webs:  codes  1,2, 6, 7, 8, 9 
|  Grain  thicknesses:  code  4 

NM 

Y 

(  Middle  webs:  code  9 
\  Slot  widths:  code  6 

NO 

M 

Outer  webs:  codes  1,7, 8, 9 

NOD 

N  ) 
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LIST  OF  SYMBOLS  (cont»d) 

Symbol 

Array 

Description 

WW 

N 

Distance  between  perforation  centers 

WW2 

N 

WW  •  WW 

GRL 

N 

Grain  length  -  U 

OD 

N 

Grain  diaaeter  -  U 

PRFD 

N 

Perforation  disaster  +  U 

PRFD2 

N 

PRFD  •  PRFD 

E 

N 

End-area  of  grain 

ALPHA  \ 

N 

Various  angles 

THETA  ) 
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